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Abstract 

We establish deep and remarkable connections among partial differential equations 
(PDEs) integrable by different methods: the inverse spectral transform method, the 
method of characteristics and the Hopf-Cole transformation. More concretely, 1) we 
show that the integrability properties (Lax pair, infinitely-many commuting symmetries, 
large classes of analytic solutions) of (2+l)-dimensional PDEs integrable by the Inverse 
Scattering Transform method (5-integrable) can be generated by the integrability prop- 
erties of the (l+l)-dimensional matrix Burgers hierarchy, integrable by the matrix Hopf- 
Cole transformation (C-integrable). 2) We show that the integrability properties i) of 
S'-integrable PDEs in (l+l)-dimensions, ii) of the multidimensional generalizations of 
the GL(M,C) self-dual Yang Mills equations, and hi) of the multidimensional Calogero 
equations can be generated by the integrability properties of a recently introduced mul- 
tidimensional matrix equation solvable by the method of characteristics. To establish 
the above links, we consider a block Frobenius matrix reduction of the relevant matrix 
fields, leading to integrable chains of matrix equations for the blocks of such a Frobenius 
matrix, followed by a systematic elimination procedure of some of these blocks. The con- 
struction of large classes of solutions of the soliton equations from solutions of the matrix 
Burgers hierarchy turns out to be intimately related to the construction of solutions in 
Sato theory. 3) We finally show that suitable generalizations of the block Frobenius ma- 
trix reduction of the matrix Burgers hierarchy generates PDEs exhibiting integrability 
properties in common with both S- and C- integrable equations. 



1 Introduction 

Integrable nonlinear partial differential equations (PDEs) can be grouped into different classes, 
depending on their method of solution. We distinguish the following three basic classes. 

1. Equations solvable by the method of characteristics PQ, hereafter called, for the sake of 
brevity, C/i-integrable, like the following matrix PDE in arbitrary dimensions [2]: 

wt + ^w Xi p®(w) + [B,w]<r(w) = 0, (1) 
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where w is a square matrix and p^\-), er(-) are scalar functions representable as positive 
power series, or like the vector equations solvable by the generalized hodograph method 

[31 HI 13 HOE]. 

2. Equations integrable by a simple change of variables, often called C-integrable [9], like 
the matrix Burgers equation [10] 

w t — Bw xx — 2Bw x w + [w, B](w x + w 2 ) = 0, (2) 

where B is any constant square matrix, linearizable by the matrix version of the Hopf-Cole 
transformation ty x = 

3. Equations integrable by less elementary methods of spectral nature, the inverse spectral 
transform (1ST) [121 H31 HH HSl US] and the dressing method pU EH [El 1201 IS] , often 
called S'-integrable [9] or soliton equations. Within this class of equations, we distinguish 
four different subclasses, depending on the nature of the associated spectral theory 

(a) Soliton equations in (l+l)-dimensions like, for instance, the Korteweg-de Vries 
(KdV) [211 12] and the Nonlinear Schrodinger (NLS) [22] equations, whose inverse 
problems are local Riemann-Hilbert (RH) problems [131 [T5] . 

(b) Their (2+l)-dimensional generalizations, like the Kadomtsev-Petviashvily (KP) [23] 
and Davey-Stewartson (DS) [23] equations, whose inverse problems are nonlocal RH 
[23 [2S] or 5 - problems [21]. 

(c) The self-dual Yang-Mills (SDYM) equation [281 EH] and its generalizations in arbi- 
trary dimensions. 

(d) Multidimensional PDEs associated with one-parameter families of commuting vector 
fields, whose novel 1ST, recently constructed in [301 [31], is characterized by nonlinear 
RH [301 EI] or 8 [32] problems. Distinguished examples are the dispersion-less KP 
equation, the heavenly equation of Plebanski [33] and the following integrable system 
of PDEs in N + 4 dimensions [30] : 

N N 

v tlZ2 - v t2Zl + ^2(v Zl ■ V £ )v Z2 - ^v Z2 - Vgv Zl = 0, (3) 

i=l i=l 

where v is an iV-dimensional vector and = (d Xl , .., d XN ). 

Each one of the above methods of solution allows one to solve a particular class of PDEs 
and is not applicable to other classes. 

Recently, several variants of the classical dressing method have been suggested, allowing 
to unify the integration algorithms for C - and S - integrable PDEs [33], for C- and Ch- 
integrable PDEs [35], and for S- and Ch- integrable PDEs [36]. In particular, the relation 
between the matrix PDE (pQ), integrable by the method of characteristics, and the GL(M,C) 
SDYM equation has been recently established in [37]. As a consequence of this result, it was 
shown that the SDYM equation admits an infinite class of lower-dimensional reductions which 
are integrable by the method of characteristics. 

In this paper we extend the results of [37], showing the existence of remarkably deep relations 
among S-, C- and Ch- integrable systems. More precisely, we do the following. 
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1. We show (in § 2) that the integrability properties (Lax pair, infinitely-many commuting 
symmetries, large classes of analytic solutions) of the C-integrable (l+l)-dimensional ma- 
trix Burgers hierarchy can be used to generate the integrability properties of S-integrable 
PDEs in (2+1) dimensions, like the N-wave, KP, and DS equations; this result is achieved 
using a block Frobenius matrix reduction of the relevant matrix field of the matrix Burgers 
hierarchy, leading to integrable chains of matrix equations for the blocks of such a Frobe- 
nius matrix, followed by a systematic elimination procedure of some of these blocks. The 
construction of large classes of solutions of the soliton equations from solutions of the ma- 
trix Burgers hierarchy turns out to be intimately related to the construction of solutions 
in Sato theory [38l [39j |40j |4J] . On the way back, starting with the Lax pair eigenfunctions 
of the derived S'-integrable systems, we show that the coefficients of their asymptotic ex- 
pansions, for large values of the spectral parameter, coincide with the elements of the 
above integrable chains, obtaining an interesting spectral meaning of such chains. It 
follows that, compiling these coefficients into the Frobenius matrix, one constructs the 
C-integrable matrix Burgers hierarchy and its solutions from the eigenfunctions of the 
S'-integrable systems. 

2. We show (in § 3) that the integrability properties of the multidimensional matrix equation 
([I]), solvable by the method of characteristics, can be used to generate the integrability 
properties of 

(a) S-integrable PDEs in (1+1) dimensions, like the N-wave, KdV, modified KdV (mKdV), 
and NLS equations (in §3.2); 

(b) S-integrable multidimensional generalizations of the GL(M, C) SDYM equations 
(in §3.3); this derivation from the simpler and basic matrix equation ([1]), allows 
one to uncover for free two important properties of such equations: a convenient 
parametrization, given in terms of the blocks of the Frobenius matrix, allowing one 
to reduce by half the number of equations, and the existence of a large class of 
solutions describing the gradient catastrophe of multidimensional waves. 

(c) S-integrable multidimensional Calogero equations j42j S3, HU HHl HE] (in §3.4). 

As before, these results are obtained considering a block Frobenius matrix reduction, 
leading to integrable chains, followed by a systematic elimination procedure of some of 
their elements. Vice-versa, such chains are satisfied by the coefficients of the asymptotic 
expansion, for large values of the spectral parameter, of the eigenfunctions of the soliton 
equations. 

3. We show (in § 4) that a proper generalization of the block Frobenius matrix reduction of 
the matrix Burgers hierarchy can be used to construct the integrability properties of non- 
linear PDEs exhibiting properties in common with both S- and C- integrable equations. 

Figure 1 below shows the diagram summarizing the connections discussed in § 2 and § 3. 

We end this introduction mentioning previous work related to our main findings, i) The 
matrix Burgers equation ([2]) with B = I, together with the block Frobenius matrix reduction 
(TT3]) . have been used in [47] to construct some explicit solutions of the linear Schrodinger 
and diffusion equations, ii) As already mentioned, once the connections illustrated in §2 are 
exploited to construct large classes of solutions of soliton equations from simpler solutions of 
the matrix Burgers hierarchy, the corresponding formalism turns out to be intimately related 
to the construction of solutions of soliton equations in Sato theory. 
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Fig. 1 The remarkable relations among PDEs integrable by the inverse spectral transform 
method, by the method of characteristics and by the Hopf-Cole transformation. 



2 Relation between C- and /S-integrability 

Usually C- and S- integrable systems are considered as completely integrable systems with 
different integrability features. In this section we show the remarkable relations between them. 

2.1 C-integrable PDEs 

It is well known that the hierarchy of C-integrable systems associated with the matrix Hopf - 
Cole transformation 

^ x = w^> (4) 

can be generated by the compatibility condition between equation (j4]) and the following hi- 
erarchies of linear commuting flows (the hierarchy generated by higher x-derivatives and its 
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replicas): 

% nm = B^d^, n,men + , (5) 

where $ and w are square matrix functions and B^ nm \ n,m G N + are constant commuting 
square matrices. The integrability conditions yield the following hierarchy of C-integrable 
equations and its replicas: 

w tnm + [w, B {nm) W {n) ] - B {nm) Wi n) = 0, (6) 

where 



(7) 



WW = W$T X) + W^-^w, n e N+, 
= I, = w, W {2) = w x + w 2 , W {3) = w xx + 2tUa.it; + ww x + w 3 , . . . . 

and I is the identity matrix. 

The first three examples, together with their commuting replicas, read: 

1. n — 1: a C-integrable iV-wave equation in (l+l)-dimensions: 

w tlm - B^w x + [w, B^)w = 0, (8) 

2. n = 2: the matrix Burgers equation: 

w t2m - B^w xx - 2B^w x w + [w, S (2m) ]K + ™ 2 ) = 0; (9) 

3. n = 3: the 3-rd order matrix Burgers equation: 



w t3m - B^w xxx - 3B^w xx w + [w, B^}(w xx + ww x + 2w x w + w 3 )- 
3B^w x (w x +w 2 ) = 0. 



(10) 



The way of generating solutions of the C- integrable PDEs (jH]) is elementary: take the 
general solution of equations ((SJ): 

r kx+ Yl BU^Hjmki A 
iir(f) = / e %?(k)dn(k), (11) 



where T is an arbitrary contour in the complex &-plane, Q(k) is an arbitrary measure and 
ty(k) is an arbitrary matrix function of the spectral parameter k, and x is the vector of all 
independent variables: x = {x,t nm ; n,m 6 N + }. Then 

iu = t^*" 1 (12) 

solves (El). 
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2.2 Block Frobenius matrix structure, integrable chains and ^-integrable 
PDEs 

It turns out that the C-integrable hierarchy of (l+l)-dimensional PDEs (jSJ), including the N- 
wave, Burgers and third order Burgers equations fl8l- ffTUl) as distinguished examples, generates 
a corresponding hierarchy of S'-integrable (2+l)-dimensional PDEs, including the celebrated 
iV-wave, DS and KP equations respectively. This is possible, due to the remarkable fact that 
eqs. and (jSJ) are compatible with the following block Frobenius matrix structure of the 
matrix function w. 



(13) 



where Im and Om are the M x M identity and zero matrices, M E N + , and j E N + 

are M x M matrix functions. This block structure of w is consistent with eqs. 01]) and (jSJ) 
(and therefore with the whole C-integrable hierarchy (jSJ)) iff matrix # is a block Wronskian 
matrix: 

^(11) #(12) #(13) 

<9^ (11) d- 1 ^ d- 1 ^ ■■■ 
d~ 2 ^ d~ 2 ^ d~ 2 ^ ■■■ 









w ® ... - 




Im 


Om 


M ■■■ 


w = 


Om 


Im 


M ■■■ 



(14) 



and 



B {im) = diag( J B (im) , B {im) ,•••), (15) 

where the blocks #( y ), i,j e N+ are M x M matrices, and B^ m \ i E N+ are constant 
commuting M x M matrices. In equations ( |T3l) -( |T5|) . the matrices w, # and B^ im ^ are chosen to 
be oo x oo square matrices containing an infinite number of finite blocks; only in dealing with 
the construction of explicit solutions, it is convenient to consider a finite number of blocks. 

Substituting the expressions ( JTBl and (fl5l) into the nonlinear PDEs (IHlTTOl. one obtains the 
following (by construction) integrable infinite chains of PDEs, for n,m E N + : 

w£l - B^wW + £ (lm) ] + [w w , B {lm) }w {n) = 0, 



w 



(n) 

t'Zra 



(16) 
(17) 



[ w V,B( 2 ™)]( w m w M + w (n) + W («+D) + [ W ( 2 ),5M] W (») + [w (n+2)^(2 m)] = 0; 

wlH - ^ (3m) (™&L + 3 ( w i^ (n) + ™t +1) ) + 3wW{wW W W + w (n+ V + «,(•*)) + (18) 
3t4 2 W"> + 3^+ 2 >) + [ W W, B (3m) ] (w£> + 2(w^w^ + w^) + W W(w<?> + 

w W w (n) + ^(n+1)) + w (2) w (n) + ^(n+2) j + ^(9 ^(n) + ^(n) + ^(n+1) j + 

[ W ( 3 ),B( 3m V (n) + M n+3) ,£ (3m) ] = 0. 

From these chains, whose spectral nature will be unveiled in § 2.4, one constructs, through 
a systematic elimination of some of the blocks the target S'-integrable PDEs. Here 
we consider the following basic examples. 
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(2+l)-dimensional TV-wave equation. Fixing n — 1 in eqs.l iToT) , and choosing m = 1,2, 
one obtains the following complete system of equations for w^\ 2 = 1,2: 

wiH - B^wl 1} + [w^,B^)w^ + [w^ 2 \B^] = 0, m = 1, 2. (19) 

Eliminating from equations (fT9l) one obtains the classical (2+l)-dimensional .S-integrable 
TV-wave equation: 

[wiH B {12) ] - [will " B^w^B^ + B^w^B^ + (20) 
[[™« B( n )]>« B( 12 )]] = 0. 

DS-type equation. Choosing n = 1,2 in equations ({TBI) , n = 1 in eq. (1171) . m = 1 in both 
equations, and simplifying the notation as follows: 

«j = tii, P = 5^ jeN +1 (21) 
one obtains the following complete system of equations for i — 1,2, 3: 

u4 1J - + [^ (2) , + = 0, (22) 

- B«u;( 2) + [tu< 3 >, B«] + [w<U, B {1) ]w^ = 0, (23) 

wff - B®w<£ - 2B^wV - 2B^w^w^ + (24) 
[wW,bW}( W WwW+wW + wW) + [w<V,bV*]w<» + [w®,B&>] = 0. 

Using eqs. f)22l) and (|23p . one can eliminate if/ 3 -* and «/ 2 ) from eq. fl24|) . In the case 
(a is a scalar), this results in the following equation for w^: 



[«,£>,£«] +a(will - B««£>B« + + (25) 

BP>w x [bV>,w] - [bW,w]w x bW) = 0. 



In the simplest case of square matrices (M = 2), with B^ = j3 diag(l, — 1) (/3 is a scalar 
constant), this equation reduces to the DS system: 

kt 2 ~ \{q„ + -p^tj - 2( V q + 2rq 2 ) = 0, (26) 
1 , 1 



Prti ~ r(r M + ^r tltl ) - 2(y?r + 2gr 2 ) = 0, 



where 



2 V /J 
- -^Ptiti + 4{rq) xx = 0, 



9 = u>2\ r = w 2i\ V = i w ii + «4 ] )«, p=—z- ( 27 ) 

a 5 



If /? = i, this system admits the reduction r = q: 



nti ~ \{(lxx + -gaStituJ - 2((^g + 2gg 2 ) = 0, (2? 
<Arx - ivtiti + 4(gg) asB = 0, 

2 _ _i 02 



becoming DS-I and DS-II if /? = — 1 and (3 =1 respectively. 
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KP. To derive the celebrated KP equation, choose M = 1, take eqs.flTTl) with n = 1,2, and 
eq. (|T8j) with n = 1, B^ = (3, B^ = — 1, where (3 is a scalar parameter, obtaining: 

«4 1} - + 2^ (1) ^ 1} + 2^i 2) ) = 0, (29) 



(1) + + 3( (wUfwW + (t0 W )a + w (D w (i) ) + 



*3 



3(^i + A« + ^Vf) + 3«;< 3 > = 0, 



where we have set m = 1 and used again the notations (I2ip . 

After eliminating and m/ 3 -*, one obtains the scalar potential KP for u = w^ 1 ', y = t 2 , 
t = t 3 : 



1 3 2 \ 3 



ut + ^u xxx + -u x j^ + —u yy = 0. (30) 
KP-I and KP-II correspond to (3 2 = — 1 and /3 2 = 1 respectively. 
2.3 Lax pairs for the 5-integrable systems 

Also the Lax pairs for the S'-integrable systems derived in §2.2 can be constructed in a similar 
way, from the system (j3J), We first observe that, due to equation (jlj), equations (0) can be 
rewritten as 

^ tnm = B^W (n) ^. (31) 

Due to the block Frobenius structure of w, it is convenient to work with the duals of equations 
(ID and (GEO): 

= -*w, (32) 
* w = _^5 (nm) W {n) . (33) 

Substituting ({TBI and (|T5|) into equations ()32H33p . one obtains a system of linear chains for the 
blocks of matrix ijf. The first few equations involving the blocks of the first row read: 

+ %(.n) w (n) + $(l(n+l)) = 0^ (34) 
£(ln) + ^(ll^lm)^) + $(l(n+l)) g(lm) = (35) 
$£0 + ^iDtfC*") + + „,(«+!)) + ^(12) j g( 2m ) w (n) + ^(l(n+2)) j g( 2 m) = q^q) 

*£? + (^ (11) 5 (3m) + 2( W W W W + w x n+ V) + wM(w<p + w (1 V n > + w^) +(37) 

w W w (n) + ^(n+2) \ + £(12)£(3m) ^(n) + ^^(n) + w (n+l)\ + ^(13)J(3m) w (n) + 
^(l(n+3)) jg(3m) _ 

where n G N + and tjKw i s the (z,j')-block of matrix 
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Lax pair for the TV- wave equation. Setting n = 1 into eqs.( l3"4"f35j) and eliminating \I/( 12 ), 
one obtains (the dual of) the Lax pair for the iV-wave equation (|20|) : 

4> tlm -4, X B^ +^ lm \w^} = 0, m = l,2, (38) 
where if) = ty( u \ The dual of it, is the well-known Lax pair of the iV-wave equation f[2"0l : 

ifj tlm - B^tp x - [£K lm ) , wW]if) = 0, m = 1, 2. (39) 
Of course, the compatibility condition of eqs. (1381) and/or eqs. ( |39i) yields the nonlinear system 

PH. 

Lax pair for DS. In this paragraph we set m = 1 in the integrable chains, and use the 
notation (I2TI) . The first equation of the dual of the Lax pair is eq. (138]) with m — 1. To derive 
the second equation, we set m = n = 1 into eq. (I36I) . and eliminate the fields \l/( 12 ), \j/( 13 ) ; using 
eq. (134"|) with n = 1, 2. In this way one obtains the dual of the Lax pair for DS-type equations: 

i> tl -^ s B^+^[B^,w^] = Q, (40) 
4> t2 + i) xx B^ + i> x [w<U, B {2) ] + $5 = 0, 

§ = [BW, W W}+W x VBW + 5(2)^(1) + [5(2)^(1)^(1). 

Therefore the Lax pair reads: 

^-sw^-taw = o, 

^ - £ (2 V- + [w®, B®)i) x - s(y)i) = 0, (41) 
a = [B^,w^]+2B^wP + [B®,wW]wM. 

The compatibility conditions of equations (140p or (141j) yield a nonlinear system equivalent to 
the system (E2H21). 



Lax pair for KP. In this paragraph we use the notations ( 1211 as well. The first equations 
of the dual of the Lax pair for KP are the scalar versions of eq. (140b ) and eq. (141b ) respectively, 
with B^> = /3 and B^ — — 1. To write the second equation of the Lax pair for KP, we must 
take the scalar version of eq. (l371) with m — n — 1, and eliminate \j>( 14 ), i = 2, 3, 4 using eq.s (l36|) 
for n — 1,2. As a result, the dual of the Lax pair reads 

1 „ 

-gK + ^xx + tyu x = 0, (42) 



and the Lax pair is 



A 3 + Ipxxx + tyxUx ~ 7^ (^J- ~ U xx ) = 0, 



^t 2 ~ ^Pxx - 2u x ip = 0, (43) 
ifka + i'xxx + 3u x if) x + ^ + u xx j ip = 0. 
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2.4 From the Lax pairs of ^S-integrable PDEs to C-integrable PDEs 

As usual in the 1ST for (2+l)-dimensional soliton equations, one introduces the spectral pa- 
rameter A into the Lax pairs (l39l) . (l4"Tl) and (j4"3l) as follows 

\xl+ J2 B^timX 

x) = x(A; x)e "• m ^ 1 , (44) 

obtaining, respectively, the following spectral systems for the new eigenfunction X : 

Xtlm - B (lm) X * - X[3 lm \ X ] ~ [B {lm) wW) X = 0, (45) 

Xtl ~ B (1) Xx ~ \[B (1) ,X) ~ [B (1) wW] X , (46) 
Xt, - B (2) Xzx + \ 2 [X, B {2) ] ~ 2\B {2) X* + [w^,B^]( Xx + \ X ) -s X = 0, (47) 
a = [BV\wW}+2BW W U + [B®,wW]wM. 



Xt2 - [3 (2X Xx + Xxx + 2 X w x ^ = 0, (48) 

Xt 3 + 3\ Xxx + 3X 2 Xx + Xxxx + 3X X w x + ^{~^ + Wxx ) = °' 

It is now easy to verify that the coefficients of the A - large expansion of the eigenfunction 
X satisfy the infinite chains (I16H18I) : 

X (A;5)=7-^— A^. (49) 

n>l 

Therefore we have obtained the spectral interpretation of such chains. In addition, since 
the infinite chains (I16H18I) for the tyW's are equivalent, via the Frobenius structure ffTSl) . to 
C-integrable systems, we have also shown how to go backward, from S- to C- integrability. 



2.5 Construction of solutions and Sato theory 

In order to construct solutions of the S-integrable PDEs generated in §2.2 from the elementary 
solution scheme (Illj) . (ll2p of the matrix Burgers hierarchy, we consider the matrices w and \l/ to 
be finite matrices consisting of n x n blocks (this can be done assuming that ^("o+i)) = 0), 
where n is an arbitrary positive integer greater than the number of blocks w^> s involved in 
the S-integrable PDE under consideration. Taking into account the structures of w and ^ given 
by eqs. (fT3|) and (fill) respectively, we have that 



w 



hi m 
0m " • • 



>M 



W 



(no) 



0v, 



M 



JM 



^-liJr(l"o) 

^-n +l^>(ln ) 



(50) 
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We remark that the no blocks j = 1, . . . , tiq of ^ are defined, via ( TIT]) , by equations 



*^')(fe)dn(A:), j = l,...,n 



(51) 



in terms of the arbitrary spectral functions while the remaining blocks are constructed 

through the equations vpw = d* -1 ^ 1 -?). Then, via (TT2T) . the components of the M x M blocks 
w^ 1 ' are expressed in terms of \l/ through the compact formula 



w 



(i) 

af3 



a(iM-M+/3) 



, a,0 = 1...,M, i 



n . 



(52) 



This formula is intimately connected to those obtained via Sato theory. To see it, we consider 
the simplest case of scalar blocks (M = 1), containing the example of the KP equation. Then 
equation f[52"j) becomes 



IV 



(0 



A 



where 



A 



9 n -2 /{ l) 9 n -2 /( 2) 



gno-1 j(n ) 



/(l) /(2) 
^o-l/(l) 0*>-l/(2) . . 

gno-j r(l) Qn -i j(2) 
Qn -i-2 j(l) Qno-i-2 y(2) . . 



./ ( 



no) 



(53) 



(54) 



Qn -i j(no) 
^n J (n ) 
Qn -i-2 j(n ) 



(55) 



j(2) y(n ) 

where f® = d~ {n °~ 1] ^\ j — 1, ... ,rio, equivalent to the formula obtained using Sato theory 



3 Relation between C7i- and S'-integrability 

Following the same strategy illustrated in §2, in this section we establish the deep relations 
between the matrix PDE (pQ), recently introduced in [2] and integrated there by the method 
of characteristics, and i) (l+l)-dimensional S-integrable soliton equations like the KdV and 
NLS equations; ii) the GL(M, C) - SDYM equation and its multidimensional generalizations; 
iii) the multidimensional Calogero systems [4"2l l4~3l |4"4"1 l4"5l H6] . In this section, matrix w must 
be diagonalizable. 

3.1 Matrix equations integrable by the method of characteristics 

Consider the following matrix eigenvalue problem 

iu(x)tf (A; x) = (A; f)A(f), (56) 
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for the matrix w(x), where A(x) is the diagonal matrix of eigenvalues, \l/(A;x) is a suitably 
normalized matrix of eigenvectors, and associate with it the following flows for \I/(A; x): 



^t mfc + E J= i^ fc P {mjfc) (A)-S( mfc )W mfe )(A) = 0, meN+, k = 1,2, (57) 

where B^ mk ^ are constant commuting matrices as in § 2 and x is the vector of all independent 
variables: x = (xjk,t nm , j,n,m e N+, k = 1,2). The compatibility between the flows ( 1571) 
implies the following commuting quasilinear PDEs for the eigenvalues: 

At mfe + £f=iA^p^(A) = 0, m eN+, k=l,2; (58) 



the additional compatibility with the eigenvalue problem ( 1561) implies the following nonlinear 
PDEs: 

N 

w tmk + J2 w x lk P {mik) {w) + [w,B^}a^ mk \w) = 0, meN + , k = 1,2, (59) 

commuting replicas of equation ([I]). 

The way of solving equations (159]) is as follows [2]. Consider the general solution of equations 
(158]) with k = 1,2, characterized by the following non-differential equations: 



A = iii/-^/) (mll) (A)t ml ,..,x m /-5]^ 1 >(A)t ml ; (60) 

\ m>l m>l 

* 12 / - J2 P (ml2) (A)^2, ■ • ■ , W - J2 P imN2) Wtm2 J , 
m>l m>l / 

where i? is an arbitrary diagonal matrix function of 2./V arguments; then the general solution 
of the linear matrix PDEs (!57j) for with the convenient parametrization ^ = 1, is given by 

= F aP (x n - P (mll) (A/3)t m i, ...,x m -J2 P (mN1 \A/3)t m i; (61) 

\ m>l m>l 

X 12 J - ]T p^iAp)^, ...,X N2 -J2 P (mJV2) (A/?)*m2 ) X 
m>l m>l 

E E Bi mfe) a(™ fc )(A^)t mfc ) 

y^fc — 1 m>l / 



> a,/? = 1,2, . . . 

where Fq,^ are arbitrary scalar functions of 2 A" arguments, with F aa = 1. Then 

w = ^Aty- 1 (62) 



solves the nonlinear PDEs (159|) . 

Now we proceed as in §2, assuming for w the same block Frobenius matrix structure (1131) . 
consistent with equations (I56p and (I57p (and then with the hierarchies (I59p ) iff the matrices 
B( mk ) are given as in (|T5|) . the diagonal matrix of eigenvalues A(x) has the block-structure 

A(x) = diag[A (1) (f), A (2) (f),---], (63) 
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and 



^r(ll) ^(22) A (2) ^(33) A (3) 2 

#(11) (AC 1 ))- 1 #( 2 2) #(33) A (3) 

#( 1 l)(A(l))-2 #(22)( A (2))-1 #(33) 

#( n )(AW)- 3 #(22)(A(2))-2 #(33)( A (3))-1 



(64) 



Using the same strategy as in § 2, we show that, 

1. if p( m -? fc ) = 0, k = 1,2, implying A = const, equations (1591) generate classical (1+1)- 
dimensional S'-integrable PDEs like the iV-wave, NLS, KdV and mKdV equations; 

2. if = (or = 0), k = 1,2, equations ([59]) generate the GL(M, C)-SDYM 
equation and its (2N + 2) - dimensional generalization; 

3. if i?( ml ) = (or a( ml ) = 0) and p(™ 2 ) = 0, equations (159]) generate Calogero systems. 

3.2 Derivation of (l+l)-dimensional 5-integrable PDEs 

Let p( mjfc ) = in (I57ti59p . implying A = const (isospectral flows), and let cr^ mfc ^(A) = A m , i.e.: 



w \$j = \]>A, 

% mk - B (mk) mk m = 0, m,keN, 



(65) 
(66) 



The compatibility conditions for the system (I65H66P yields, for m € N + : 

w tmk + [w, B^]w m = 0. 



(67) 



We remark that these equations are equivalent to eqs. (JHJ) with d J x w = 0, Vj. Consequently, the 
discrete chains generated by the eq.( |67j) . with m = 1, k = 1, 2 and m — 2, 3, k — 1, are given 
by the eqs. flUDHHJ with d j x w = 0, Vj: 



w 



H + f^ 1 ),^)] + [w^,B^]w^ = 0, fc = 1,2, 



(68) 



<° + B< 2 )] (icWtflW + u/"^) + [«;( 2 ), S( 2 V (n) + [^ (n+2) , 5 (2) ] = 0, (69) 



(n) 



t 3 
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+ B( 3 '] W W( W (1, w (n) + w (n+1) ) + w (2) » (n) + w (n+2) + 



(70) 



^\m\(w^w^+w^ n+1 A + [w^,B^]w^ + [w^ n+3 \B^} = 0, 

for n G N+, where, in equations (!69l . (170]) . we have used the simplifying notation ([2"T]) . 

(l+l)-dimensional A^-wave equation. Setting n = 1 in eqs. (168]) . and eliminating w^ 2 \ 
one obtains the well-known S- integrable A^-wave system in (l+l)-dimensions: 



(71) 
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NLS. Eg. ( 1681) with k = 1, n = 1,2 and eq. (l69i) with n — 1 are a complete system of PDEs 
for w®, j = 1,2,3: 

w$ + [ W W,BW] + [ w V>,bW]w<V = 0, (72) 
u£ 2) + [w®, £«] + [ w M,BW]wW = 0, 

wg> + [wW,BM](wU W W+wW) + [wV>,B®]wM + [w^,B^\ = 

In the case B® = aB<» (a is a scalar constant) this system results in the following equation 
for w^ 1 ': 

[wg\BW] + wl% - «K ( > (1) ,5«] +a{[w^,B^]w^) t = (73) 

If, in addition, 

M = 2, = diag(l,-l), (74) 

this equation yields the celebrated NLS system: 

1 i 

~Qt 2 ~ -q-nn ~ 4? V = 0, (75) 
a 2 

~ r t 2 + 7T r nn + Mr =0 
a 2 

for the off-diagonal elements of : q = , r = if% ■ The NLS equation iq t2 + \q T1T1 +4q 2 q = 
corresponds to the reduction r = q, a = i. 

KdV and mKdV. Eq. (l68l) with k = 1, n = 1, 2, 3, and eq. (17Dl) with n = 1 yield a complete 
system of PDEs for j = 1, 2, 3, 4: 

4 1} + [u,< 2 >, £W] + [tcW.BWjujW = 0, (76) 
+ + [«;« SW]^ 2 ) = 0, 

^ + [w« ^(^(^V 1 ) W 2) ) + tu< 2) «> W W 3) ) + 

^(^W 1 ) W 2) ) + [w^^V^ + [™< 4 >,b< 3 >] = 0. 

In the case B^ = —B^\ this system reduces to the two coupled matrix equations 

+ [u;^ + [wW,BW]wW = 0, 

= wi% + [5« ) 4 1 V (1 ^ (1) +M) (2) )+4 2 ^ (1) ] - ([B^W^ . (77) 

If, in addition, the choice (1T4"|) is made, the system (1771) becomes the mKdV system: 

% + ^ftititi + = 0> ( 78 ) 

1 

r * 3 + | r *ititi + 6r tl rg = 0, 

where q = w^, r = «4i\ reducing to the KdV equation g t3 + l?*^^ + 6gg tl = and to the 
mKdV equation % + -xq^tx + 6g 2( Zti = if r = 1 and r = q respectively. 
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3.2.1 Lax pairs for S'-integrable PDEs in (l+l)-dimensions 

As in §2, in order to derive the Lax pairs for the above S'-integrable PDEs in (1 + 1) -dimensions, 
it is convenient to write the system (I65f66l) in the equivalent form 



and consider the dual system 



w \& = \]>A, 

^ _ B {mk) w m V = 0, m e N H 



# t , + ^B {mk) w {m) = 0, m E N H 



(79) 
(80) 



(81) 
(82) 



Taking into account the block structure of the matrix w given by eq. fflUj) and considering the 
first rows of eqs. (l81|l82l) . we obtain the following spectral chains, for n 6 N + : 

Vjr(ln)£ = ty(U) w {n) + $(l(n+l)) ^ (gg) 
Vjy(ln) + §(11) g(lt) w (n) + §(l(n+l)) £(lfc) = q (84) 

^lfc * ' 

§(ln) + £(ll)£(2fc)( w (l) w (n) + ^(n+1)) + §(12) £(2*)^) + §(l(n+2)) 5(2*) = 0> ( g 5) 
§gn) + $(11)£(3*) ^(D( w (\(«) + w (n+l)) + + w (n+2)^ (86) 

+ §(12) £(3fc) + y;(n+l)^ + §(13) B {Zk) W {n) + §( 1 ( n+3 )) 5< 3fc ) = 0. 

where £ = A^. Setting n = 1 into eqs. (l83|l84|) and eliminating §( 12 ), one gets the dual of the 
Lax pair for the (l + l)-dimensional A^-wave equation f)7ip (tp = tyt 11 '): 

$ hk + BW lh) + i>[B {lk) ,w {1) ] = 0, k = 1,2. 



57) 



Eq. (1871) with k = 1, written in terms of (1211) . is the first equation of the dual Lax pair also for 
eqs.([72D and (|76|). 

Setting k = n = 1 into eq. fl85|) and eliminating \|>( 12 ), vjH 13 ) using eq. fl83|) . one gets the second 
equation of the dual Lax pair for (1721) : 

(88) 



The second equation of the dual Lax pair for the eq. (176|) results from the eq. (1861) . k = n = 1 
after eliminating §w), j = 2,3,4 using eq. (!83|) . In view of conditions (1741) the complete dual 
spectral system reads: 



4 - *fy 



-2qr 



1 
-1 

1 
-1 

2qr 



2q 
-2r 

2q 
-2r 



0, 

+ 



(89) 



4g 2 r 



4r 2 g gr T1 — rg T1 



J_ y /I /k>2 



2 ' nn 
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The duals of the equations ( 1H7H891) read: 



ifj hk -3 lk ^S- [B (lk \w (1) }ij = 0, k — 1,2, (90) 
i) t2 - B^^S 2 - [B {2) ,w (1) }^£ -sif; = 0, (91) 



Ai + 



A 3 



1 o 

-1 

1 
-1 



-2qr -q Tl 
—r T1 2qr 



ip£ - 
tp£ + 




-2r 



-2r 



2q 


2q 




V = o, 



(92) 



i 



rq Tl - qr Tl ~g T1T1 + 4g 2 r 
-hr-nn - 4r 2 q qr Tl - rq Tl 



if} = 



Eqs. (!90|) are the Lax pair of the A- wave eg. (171 p . eq. (!90j) with A; = 1 and eg. (191]) are the 
Lax pair of eq. (j73l) . reducing to the NLS system if ( 1741) holds, and eqs.( l92l) are Lax pair of the 
system (1781) reducing to either KdV or mKdV. 



3.2.2 Prom the Lax pairs of (l+l)-dimensional S'-integrable PDEs to the Ch- 
integrable eqs. (1671) 

We show that the eqs. (1671) may be derived from the spectral problems obtained in §3.2.11 Let 

J2 (-l) i B l - ik h ik A i 

ip(A;x) = x(A;x)e 1 '^ 1 . (93) 

Then equation pH2]) yield: 

Xt lk +A[^ lk \x}-[B (lk \w^]x = 0, A; = 1,2, (94) 
Xt2 - A 2 iB {2) , X] - [w^,B^]A X -s X = 0, (95) 
s= [B<- 2 \w®] + [BW,wW]wW. 



Xti + 
Xt 3 - 





" 1 












-1 


,x 




" 1 












-1 


,x 



A- 

A 3 + 



2q 
-2r 

2q 
-2r 



-2qr —q 
2qr 



n 



n 



XA + 



rq 



n 
2 r m 



X = 0, 
XA 2 + 



(96) 
(97) 



qr Tl % q Tin + 4q 2 r 
— 4r 2 q qr T 



n 



X 



It is now easy to verify, as in §2, that the coefficients of the A - large expansion of the eigen- 
function x satisfy the infinite chains ( 1681) . ( 1681) . (1691) . ( T70l : 



X (A;y)^/-^(-l)^A^, 
i>i 



(98) 



clarifying their spectral meaning. In addition, recompiling such coefficients into the block Frobe- 
nius matrix, we reconstruct the matrix equations (1671) from the Lax pairs of (l+l)-dimensional 
S'-integrable PDEs. 
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3.3 Derivation of the SDYM equation and of its multi-dimensional 
generalizations 

Now we take, in equations (J56I) , (J57J) , = (or = 0), m = 1, p^ k \A) = A*, t n = t, 

= r, Xji = Xj, Xj2 = jjj, obtaining the system 

wty = (99) 

N 

3=1 

N 
3=1 

whose compatibility condition yields 

TV 

wt + ^w^w* = 0, (100) 

3=1 
N 

W T + '^2Wy j W j = 0. 

3=1 

We proceed as in the previous sections but, before considering the derivation in the general 
case, quite complicated, we illustrate the simplest two examples. 

3.3.1 N = 1: the GL(M, C) - SDYM equation 
The compatibility condition of the system fl99l) yields 

w t + w Xl w = 0, w T + w yi w = 0. (101) 
Let w and ^> be given by eqs. ({TBI and (1641) respectively. The first rows of the matrix equations 



generate the chains, for n G N + [37] : 

wl n) + w^w {n) + = 0, (102) 



wi n) + wftw™ + w^ +1) = 0. 



Setting n—1 and eliminating w^ 2 \ we derive the well-known GL(M, C) — SDYM equation: 

(1) -«•«+[<>,«;«] = 0. (103) 



To derive the Lax pair of (11031) . we write first the dual of system (1991) in the convenient 
form: 

= A*, (104) 
* t + tf^w = 0, # T + ^ yi w = 0. 

Using again (Tl3|) and (I64p . the first rows of equations (j!04p appear in the form: 

#( n V*) + = ^^(^), (105) 

^ ln) + ¥£) W W + ^(l(n+l)) = 0> 



lf n > + ¥}^w^ + ¥}} n+1 » = o, 



yi ' yi 
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where £ = Setting n = 1 in eqs. (11051 ) and eliminating \|>( 12 ), one obtains the dual of the 
Lax pair of fll03j) for the spectral function tp = xp^ 11 ); 



$ t + (e$) xi = (106) 



Then the Lax pair of ( 11031) reads: 



iPt + ip xi £ + w^iP = 0, (107) 
Vv + 4j Vi £ + u>£ty = 0. 

Vice-versa, it is easy to verify that the coefficients of the £ large expansion of the eigen- 
f unction -0 in ( 110 7ft are the elements of the chain ( 11021) : 

-> / - (108) 

i>i 

obtaining the spectral meaning of such chains. As a consequence, one reconstructs eqs. dlOip 
from the Lax pair (11071) of the SDYM equation. 

3.3.2 N = 2: a generalization of the SDYM equation in 6 dimensions 

The compatibility condition of the system (1991) yields now 

w t + w Xl w + w X2 w 2 = 0, w T + w yi w + w y2 w 2 = 0. (109) 
The corresponding chains read: 

wl n) + wjgwM + + + «/ n+1) ) + «4 2 > (n) + 4" +2) = °> ( n °) 

W W + w$wM + + wg\w^w^ + w^) + w®wW + w£W = 0. 

Setting n — 1 and n = 2 in (IllOp and eliminating the fields w^ 3 \ one obtains the following 
integrable system of two nonlinear PDEs in 6 dimensions for the fields w^\w^: 

(1) (1) , (2) (2) . (1) (1) (1) (1) . (1) (1) (1) (1) . r (2) (1), 

W X2 ' T - Wy 2 ' t + W X2 ' yi - W Xl y 2 + W X2yi W (L> - W Xl y 2 W (L> + W yi 'W X2 - W Xl 'Wy 2 + [Wy 2 , W X2 \ 



V / W2 V / j/2 



w { St - w { y^ t + [iu£\ wi 1 }) + w x % - w y 2 2 ] + [w$,w2?] + [wyW w { X2 ] ] + [w { X2 \ w { y^}w {1)2 + 
(w x % 2 w^ - w x lM l) + w x % 2 - wM + u,W (wP _ W W W M + [„£>,„,«]) _ 

(111) 

reducing to ( 11031) for w^ 1 ' if the fields do not depend on £2,2/2- 

To derive the Lax pair of ( 11111) . we consider again the dual of system ( 1991) : 



^w = A^, (112) 



* t + V xi w + ^ X2 w 2 = 0, ^ T + + ^ 2 w 2 
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Using ( fl3l) and ( |64"1) . the first rows of equations ( 1112ft appear in the form: 

+ = (113) 

§J ln) + ^)«,W + + ^(a/WiflW + w^ 1 )) + + *£ (n+2)) = 0(114) 
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(n+1) 



) + *£ a M»> + m n+2 » = 0(115) 



Setting n = 1 in equations (j!14|115p . and eliminating \j>( 12 ) ; \]/( 13 ) using eqs. (11131) for n— 1,2, 
one obtains the dual of the Lax pair of (11111) for the spectral function ip = vjK 11 ); 



a) 
i 2 ' 



(116) 



3=1 

2 



3=1 

Then the Lax pair of (11111) reads: 

ijH + i> Xl £ + iP X2 £ 2 + (w% + wgw® + w£) i> + wjgipe = 0, (117) 
Vv + i> yi E + ^ y2 £ 2 + (w$ + wftw® + w£>) i> + w^iPS = 0. 

As before, it is easy to verify that equation (11081) holds, namely that the coefficients of the 
£ large expansion of ip in (11171) are the elements of the chain (IllOp . Therefore one reconstructs 
equations (11091) from the Lax pair (11111) of the six dimensional generalization (II lip of the 
SDYM equation. 



3.3.3 Multidimensional generalization of the SDYM equation 

Motivated by the above formulae for the simplest cases N = 1,2, here we discuss the general N 
situation. If w is the block Frobenius matrix (fTBl) . the power exhibits the following structure 

yjirM y)(i;i3) 



w' 



yjUv 1 ) u)0';i 2 ) yjird 3 ) 
Im 0m 
®m Im 0m 



where the matrix blocks w are defined by the equations 



3-1 



w 



(n)x 



1=1 



2<k<j, 

in terms of the basic blocks j > 1. The first few examples read: 



w {j ' jn) = w {ny 



^(3;ln) = (tt (l) )2w (n) + + ^(n) + 

^(4;ln) = w (l) ^ w (D)2 w (n) + w (D w (n+l) + ^(^(n) + ^(n+2) \ + 



;ii8) 



(119) 
(120) 

(121) 
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Furthermore, evaluating the (ln)-block of u> J+1 , written as (w^w), we obtain the additional 
formula 

^0-l;l(n+l)) = -(,-;!„) _ l5 (M;il) w («) ] j > ^ „ > ^ ( 122 ) 

implying 



n-l 



_ ^ T2 (j+»-i;ii) tl ,(n-0 > j,n > 1. (123) 



Eq. (|123p reduces, for j = 1, to the following equation: 

j'-i 



i") = ^ + j > 1 (124) 



i=i 



useful later on (compare it with equation (I119p for n — 1). 

Using equation (ITBl . the system (llOOp generates the discrete chains 



w 



N j 

i n) + E (» : '"" '• |:: " : + 4f ->) = °> ( 125 ) 



i=i i=i 

iV , j 

j = l ' i=l 



Setting n = 1, . . . , JV in equations (11251) . one obtains a determined system of 2N equations for 
the fields w {i) , % = 1, . . . , 2N. As in the previous two illustrative examples for N = 1,2, it is 
possible to eliminate the N fields i = N + l, . . . , 2N, obtaining a system on N equations in 
(2iV+2) dimensions for the remaining fields i = 1, . . . , N. Such system, which provides the 
natural multidimensional generalization of the SDYM equation fl 1 3 j) . is conveniently written 
as follows 

p (N,0)_ q (N,0) + [q{0)jp{0)]=0i 



q (N,n) + £ / (N,n-j) _ q (N,n-j)\ + £ ^N,) ^N,n-^ = ^ j < „ < ^ _ ^ (126) 

i=i v y i=o 

where the fields p( N >i\ q\ N ^> are suitable combinations of the iV fields w^ n \ n = 1, . . . , N: 

N /s—j—1 \ 

«=i+i V l=i J 



s=j+l \ 1 = 1 
N fs-j-1 



8=j+l \ 1=1 



(127) 



the first few read as follows: 

JN,N-1) = W (X> 



p {N,N- 


1) 


— UJx N , 




p (N,N- 


-2) 


= WxLi 


+ wfl 4 


p {N,N- 


-3) 


-w {1) 


+ will 


q (N,N- 


•3) 


-w (1) 


' w yN-i 



Vn-i 



;i28) 



45 + (41. + 45) + 45 (> )2 + ^ (2) 

To show it, it is more convenient to go through the Lax pair derivation. 
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Lax pair. The system dual of ( 1991) reads 

= Af , (129) 

TV TV 

3=1 3=1 
and it is conveniently rewritten in the equivalent form 

^Iw = A*, (130) 

w y y- s ~ x = 0(132) 

As before, the block (1, 1) of the matrix equations (1131l) .( fT32l leads to the dual of the Lax pair 
{£ = A^) of the multidimensional SDYM equations: 

TV TV-1 

j> t + ^2(S^) Xj = J2 (133) 

3=1 j=0 
TV N-l 

fa + = E &W NJ) , (134) 



TV 

* t + E( Ai *k 


JV 

-E*(™% 


TV 

= ¥ t + E( A ^k 


TV-1 

- ^A a tf 


TV 

E 






i=i 




3=s+l 


TV 


TV 


TV 


TV-1 


TV 




3=1 


3=1 


- ^ A S ^ 

s=0 


E 

j=8+l 



3=1 j=0 

where p( N >fi and gWi) are defined in terms of w;W and their derivatives in (j!27p . 
Then one derives the corresponding Lax pair 

TV TV-1 

^ + E + E ^ V^' = 0, (135) 

3=1 3=0 
TV TV-1 

+ E + E q (N ' j) i>£ j = o, (136) 

3=1 3=0 

together with its compatibility condition, the following determined system of 2iV equation in 
(2tV + 2) variables for the fields p iN ' j \ q (N > j \ j — 0, . . . , N — 1: 

p(iv 1 o)_ g (iv ) o) + [g (o) )p (o) ] = 0! (137) 

n n 

pf •») - ^ + £ (pg^ - gW-4>) + £ [g^.pW-^] = 0, 1 < n < N - ^38) 

3=1 3=0 
TV TV-1 

E ^ -'/.:"•" ") • E [? (JVj ' ) ,P Wn - J ' ) ] =0, N<n<2N-2, (139) 

j=n-TV+l j=n-TV+l 

We remark that only the first A" equations (j!37p . (1138p involve derivatives with respect to the 
time variables t,r; the remaining N equations (jl39p . (ll40j) . providing a set of relations among 
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the 2N fields p^ N ^\q^ N ^\ are remarkably parametrized by equations (11271) in terms of the N 
fields j = 1, . . . , N. Therefore one is left with equations (I126p . (ll27p . 

We also remark that the generalization (I137H140]) of the SDYM equation is known in the 
literature [HI 08] , to be generated by the Lax pair 



N N-l 



ifit + Y, A ^ + 2 Xp im il> = 0, (141) 

j=l j=0 
N N-l 

ipr + J2 X ^Vi + H A^W'V = 0, (142) 



differing from (I135|ll36p by the fact that here A is just a constant and scalar spectral 
parameter. 

Therefore the remarkable derivation of (1137111401) from the matrix equations (11001) and its 
integration scheme has allowed one to uncover the following two important properties of the 
system fil37til40l) . 

• Half of the equations of the system (1137111401) (the non evolutionary part) can be parametrized 
in terms of the blocks j — 1, .., N of the Frobenius matrix w, reducing by half the 
number of equations. 

• Equations (I137H1401) turn out to be associated with the novel Lax pair (I135II136I) . in which 
the diagonal matrix £ satisfies, from (158|) . the integrable quasi-linear equations 



N N 

St+T, = °, £ r + E V = 0. (143) 

Therefore, as it was already observed in [37] in the case of the SDYM equation (I103p . the 
integration scheme associated with such a novel Lax pair makes clear the existence of a 
rich solution space exhibiting interesting phenomena of multidimensional wave breaking. 
A detailed study of these solutions is postponed to a subsequent paper, together with 
the comparison with the finite gap solutions of the SDYM equation constructed in [49] , 
and associated with a Riemann surface with branch points satisfying equations (I143P for 
N = 1. 



From the Lax pair of the multidimensional SDYM to the integrable chains ( 11251) 

As for the particular cases iV = 1, 2, in this section we show that the £-large limit of the Lax 
pair (I135II136I) yields the expansion (11081) for the eigenfunction ip. Therefore the coefficients of 
the £-large expansion of the spectral function associated with the S'-integrable multidimensional 
generalization of the SDYM equations are solutions of the nonlinear chains (11251) . providing the 
spectral meaning to such nonlinear chains. In addition, recompiling the matrices into the 
block Frobenius matrix w, via (TT3"|) . one establishes a remarkable relation between the Lax pair 
of the multidimensional SDYM and the basic matrix equations (11001) . solvable by the method 
of characteristics. 

To show the validity of the expansion (11081) . we substitute it into the Lax pair (I135|136l) . 
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obtaining the following pairs of equations 

N-l 

I AT i~\ 

P 



(m = j2 ; ' + P ^) w u-A + v £w i o<i<N-i, 

j=i+l ^ ' 

W) = J2 (w^ + qWwU-Q) + w [ y N N - l \ < i < N - 1, 
j=i+i v / 

4° + ^ (<4 + ° + P (JVJ) ^ (i+i) ) + wi N N +l) + pW°>«;(0 =0, i > 1, 
i=i V 

iV-1 



(*) , 



V-l ' ' ( 145 ) 

£ + q^w^) + w ( y N N +i) + qWw® = 0, % > 1, 

corresponding respectively to the condition that the coefficients of the positive and negative 
powers of £ are zero in all orders. Equations (I144p are identically satisfied using the definitions 
(11271) and (11191) of p( N >i\q( N >j> and w^' u \ To show that also equations fll45j) are identically 
satisfied, we compare them with equations (11251) and use again (11271) . (11191) . to finally derive 
the equation 

N k / k—l k—l — 1 \ 

S 2 W S ( 2 «c^ (w - i+liln) - J] * (fc -^' ;11) ^ 0+ ") )=0, n > 1 (146) 

fc=l /=1 \i=l 3=0 / 

identically satisfied, due to ffTT9|) and ffl~24"j) . 



3.4 Derivation of Calogero systems 

In the Lax pair of the S-integrable Calogero systems [121 021 HH SSI HE], the spectral problem 
is 1-dimensional (like that of, say, KdV and NLS), while the equation describing the evolution 
of the eigenfunction is a multidimensional PDE in which an arbitrary number of additional 
independent variables are graded by powers of the spectral parameter; in addition, such spectral 
parameter satisfies a quasilinear PDE. 

It is therefore clear that Calogero systems combine properties of the (S-integrable PDEs in 
(l+l)-dimensions of §3.2l with properties of the multidimensional generalizations of the SDYM 
equation of §3.3 and, to generate them from equations (|59|) . we have to start with the matrix 
equation (156]) and with evolutionary system (|57|) in which p^ 1 ' = cr( m2 ) = 0: 

= A, (147) 

% ml -B^ ml ^a^(A) = 0, meN + , 

*u + 5ji*^(A) = 0- ( } 

Here we illustrate the construction of the simplest example of Calogero system, corresponding 
to N = m = 1 and a^ u \A) = p< 112 )(A) = A, B = using the notation B = r = t n , 

t = t 12 , x = X\2- Then the compatibility between the two equations (I148j) yields the following 
quasi-linear PDEs 

A t + A X A = 0, A r = (149) 

for the matrix of eigenvalues; the compatibility between equations (I147P and f!148j) yields instead 
the following matrix equations for field w 

w T + [w, B}w = 0, (150) 
w t + w x w = 0. 
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Using the Frobenius structure ( 1131) of w, we get the following discrete chains, for n e N+: 

w (n) + [ ty ("+i) ) b] + [^W, %W = 0, (151) 
wi n) + wPw^ + W ^+V = 0, 

coinciding with eqs. (l68l) for k — 1, and with (1102a ). Fixing n — 1 in (11511) . and eliminating 
i// 2 -*, one gets the following matrix PDE 

K {1) , B] + [ujWioW, S] = w£ + (V 1} , . (152) 

Let M = 2 and 5 = /?diag(l, -1); then eq. (TT52]) yields 

-PQt ~ \q,r - ^ 2 qd;\qr) x = 0, (153) 
Pr t -^r„-4p 2 rd; 1 (qr) x = 0. 

If r = q, /3 = i, then the above system becomes the following (2+l)-dimensional integrable 
variant of NLS: 

iqt + ^Qxx - Md^iqq)* = 0, (154) 

studied in [50l [51] . 

Using the dual of equations (!82|) with m = k = 1, and equations (I104a .b). we derive the 
dual Lax pair ((E3) for k = 1 and (fl06k )) for (fl52l) : 

^ T + ^ J B + ^[ J B,w (1) ] = 0, (155) 
& + = ^ 

and the corresponding Lax pair ( (190]) for A; = 1 and (1107a )): 

ip T - Btp£ - [B, w (1) ]ip = 0, (156) 

V>* + + = o 

for the system (11521) . 

We end this section remarking that integrable PDEs associated with Lax pairs with varying 
spectral parameter have been studied also elsewhere, see, for instance, [52J. 

3.5 Construction of solutions 

The construction of solutions for the three classes of S'-integrable PDEs derived in §3.2, 3.3 
and 3.4 from eq. (l59|) . is based on the solution of the algebraic equations (I60II61I) . taking into 
account the block-matrix structure (1641) of Then the independent blocks ty[ u \ i = 1,2, . . . 
are characterized by the explicit formula: 



iT/(") _ rp( u ) 



(x n -J2p {mll) (4 ) )tmi,...,x m - Xy^HA^tna; (157) 

\ m>l m>l 

X12 - X P (ml2) (A?)t m2 , ...,X N2 -J2 P {mN2) ^f)tm^\ X 
m>l m>l 

E E B { ™ k ^^){\f )t mk 



, a. 3 — 1. .... M. i-\,2,... 
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where F^ % ) are arbitrary scalar functions of 2N arguments (such that Faa =1)? while the re- 
maining blocks are given by the equations = *(#)(A0'))i-\ nee A and ^ are constructed 
in this way, the blocks are obtained, from eq. fl62|) . through the compact formula: 

^^(^A*- 1 )^.^), a,(3=l...,M. (158) 

In the case of the (l+l)-dimensional S'-integrable models of § 13.21 corresponding to p(^ fc ) = 0, 
A is an arbitrary constant diagonal matrix and formula (11571) reduces to: 

*S = *S ) e* =lm - 1 > a,P=l,...,M, i = l,2,..., (159) 

where are constant amplitudes. Then the solution (11581) is a rational combination of 
exponentials. It would be interesting to compare, in this case, the solution space generated by 
(11581) with that generated by Sato theory. 

In the case of the Calogero systems of § 13.41 when 

p (mjl) = a (m2) 

= 0, formula (11571) reduces 

to: 



P (ml2) (A?)t m2 , ...,x m -J2 P (miV2) (A?)t m2 ) x (160) 

m>l m>l / 

E Bi ml V( ml )(A«)i ml 
e"^ 1 , a,/3 = 1,...,M, i= 1,2,.... 

where F^g are now arbitrary functions of N arguments and A is the implicit solution of the 
nondifferential equation 

A = i?(x 12 /-^p( ml2 )(A)t m2 ,...,W-^P (miV2) (A)t m2 ) , (161) 

\ m>l m>l J 

following from the eq. (l60]) . where E is arbitrary diagonal matrix function of N arguments. 



4 Generalizations 



The block Frobenius matrix (|T3j) is not the only possible structure of w allowing one to generate 
new types of integrable nonlinear PDEs, starting with C-integrable and C/i-integrable equations 
([6]) and (1591) . A more general representation is given by the following block matrix: 

" • • • " 

W (21) W (22) . . . 



W 



(162) 



where each block has one of the following two block matrix forms: 









Im 


Om 


Om 


Om 


Im 


Om 




w (»i;2) 


y;(y;3) 


Om 


hi 


Om 


Om 


Om 


hi 



(163) 
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and the blocks w^' k ' are M x M matrices. 

To provide consistency of this structure with the Lax pairs (I4"f5]) or (I56ll57p . we must take 
appropriate structure of matrices £>( nm ) and Consider the simplest example of eqs. (|4|5|) with 
n = 1. In this case w satisfies the N-wave equations ([8]). Let 



w 



5,(11) C (12) 
£1(21) C (22) 



(164) 



Then I?( lm ) must be taken in the form 

£( lm ) = diag(-B (ml) , _B (m2) ), B {mj) = diag{B {rnj) , B {mj) ... ), j = 1, 2, (165) 
where B^ are diagonal matrices, and \1/ must have the following block structure: 



#(21) ^(22) 



(166) 



where vpw) must satisfy the following system of linear PDEs (consequence of eq.fllj)), for i > 2: 

= ^(H;(i-l)i) + ^(21;^)^ ( 167 ) 
^(ia;«) = ^(12i(t-l)j) + ^(22;ij) } (168) 
^(21;ij) _ ^(ll;y) + ^(21;y)^ ^gg) 
^(22;ij) = ^(12;y) + ^(22;ij)_ 

In view of (11641) . eq.® reads (£j = tij): 

5 (11) _ 5(^)5(11) + ^(llJ^Wjfi'dl) + (C7(12)s(i2) _ 5(il) C (12) )C r(21) = 0> (in) 
c fl2) _ ^(0)^(12) + ^(11)^(^1)^(12) + ^(12)^2) _ ^(,1)^(12)^(22) = ^ 
C[ 21) - i^2)(7(21) + ^(21) 5(a) _ g{i2) C {21) ^(11) _j_ ^(22)^(12)^(21) = 
£,(22) _ £(i2) C (22) + (£(21) _ £(i2) £,(21)^(12) + ^(22) ^(i2) j ^(22) _ 

where z = 1,2. Writing the block (11) of each of these equations, one obtains 

- B^w^ + [wW, BW] W W + [w^ 2 \BW] + (u/ 12 ; 1 )^ 2 ) - (172) 
BWw^w&M = 0, 

_ bW w £W + [u,^ 1 ),.^ 1 )]^ 1 * 1 ) + (w (12;1) 5 (j2) - B^wW))^ 1 ) = 0, 

w (21;l) _ £02)^(2151) + ( w (21;l)5«l) _ £(0)^(21;!) )w (ll;l) + ^ (21; 2)£(il) _ ^^l;^ + 
[ w W) )j g(i2)] u; (21 ! l) =0j 

w (22;l) _ ^(0^(22;!) + ( w (21;l)£(tt) _ fra)^!;!)^!) + ^1)^(0)^(22-1) = Q 

Eliminating k/ 11;2 ) and it/ 21 ' 2 ) from the eqs. (1172k . c) with i = 1,2 and taking eqs.( j!72b .d) 
with % = 2, one obtains the following (2+l)-dimensional evolutionary system of PDEs in the 
time variable t 2 : 

[£(1)^(21)1 _ [£(2) ? 5(11)1 _ q 5 (173) 
q t2 - B^q x + [v, B^]q + (qB^ - B^q)u = 0, (174) 

£(1)5(21) _ £(22) £(1) _ 5(2)3(11) + 5(12)^(2) = Q) (175) 

u t2 - B^u x + (wB^ - B™w)q + [u, B {22) }u = 0, (176) 
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where 



and 



£« = v u - B( fl )» ( + [v, B( a )]<u + (gi^ 2 ) - B^q)w, 2 = 1,2, 
£(*) = - B^w x + (wBW - B (l2) w)v + [m, B( i2 )]to, i = 1, 2, 



(177) 



u = «,(i 1 ; 1 ), g = «,( 12 ^, «, = «,( 21 ;D, u = w^ 1 \ (178) 

supplemented by eqs. (1172b . d) with i — 1: 

% - + [v, B (11) ]g + (g5( 12 ) - B^q)u = 0, (179) 

M tl - B (12) M;E + («;B( n ) - B( 12 )«;)g + [u, B^]u = 0, (180) 

that can be considered as compatible constraints for the evolutionary system (I173p - (I176I) . 

The Lax pair, as well as the solution space for this system, can be obtained following the 
procedures described in §2.31 and §2.51 

If w = q = u = 0, one obtains the S'-integrable (2+l)-dimensional A^-wave equation fl 1 73[) 
for v (see eq.f l20l) ). Instead, if v = w = q = 0, one obtains the C-integrable (l+l)-dimensional 
A^-wave type equation (11761) or (I180p for u (see eq.flS])). Therefore equations (jl73j) - (jl76l) and 
(I179lll80p can be viewed as nonlinear PDEs sharing properties of S- and C- integrable systems. 
One can show that this property is shared by all the PDEs generated by reductions of the type 
(I162II163P ; recovering, in particular, the (n+l)-dimensional (n > 2) nonlinear PDEs constructed 
in [M] by a generalization of the dressing method. 



5 Summary and future perspectives 

We have established deep and remarkable connections among PDEs integrable by the inverse 
spectral transform method, the method of characteristics and the Hopf-Cole transformation. 
These relations can be used effectively to construct, for the generated S-integrable PDEs, the 
associated compatible systems of linear operators, their commuting flows and large classes of 
solutions. These results open several research perspectives. 

1. Use of the above derivation of the S - integrable systems to investigate the corresponding 
space of analytic solutions generated from the seed solutions of the original C - and Ch - 
integrable PDEs. In particular, 

(a) the connections between such solution space and that generated by Sato theory. In 
the KP case, the two solution spaces coincide; in other cases the connection is, at 
the moment, less clear. 

(b) The use of the quasi-linear PDEs for the eigenvalues to study in detail the wave 
breaking phenomena associated with solutions of the SDYM equation and of its 
multidimensional generalizations. 

2. Search for the integrable systems that should generate, through a suitable matrix reduc- 
tion, the integrable PDEs equivalent to the commutation of vector fields, like equation 
([3]), a class of S'-integrable systems not fitting yet into the general picture illustrated in 
this paper. 
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3. Generalization of the techniques presented in this paper to generate novel integrable 
systems, possibly in multidimensions. In particular, a systematic use of group theory 
tools to explore reductions different from the block Frobenius matrix one. 

4. Construction of the discrete analogue of the results of this paper. In this respect, we 
remark that, while the discretization of the results of §2 does not present, in principle, 
any conceptual problem, and will be the subject of a subsequent paper, the discretization 
of the results of §3, if p^ k ^> ^ 0, is not clear, since a satisfactory discretization of the 
method of characteristics and of equation ([1]) for pW ^ o, in the scalar and matrix cases, 
are, at the moment, unknown. 
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